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Abstract. Cooling simulations of neutron stars (NSs) are performed assuming that stellar cores consist of neutrons, 
protons and electrons and using realistic density profiles of superfluid critical temperatures T cn (p) and T cp (p) of 
neutrons and protons. Taking a suitable profile of T cp (p) with maximum ~ 5 x 10 9 K one can obtain smooth 
transition from slow to rapid cooling with increasing stellar mass. Adopting the same profile one can explain 
the majority of observations of thermal emission from isolated middle-aged NSs by cooling of NSs with different 
masses either with no neutron superfluidity in the cores or with a weak superfluidity, T cn < 10 8 K. The required 
masses range from ~ 1.2 Mq for (young and hot) RX J0822-43 and (old and warm) PSR 1055-52 and RX J1856- 
3754 to ~ 1.45 Mq for the (rather cold) Geminga and Vela pulsars. Observations constrain the T cn (p) and T cp (p) 
profiles with respect to the threshold density of direct Urea process and maximum central density of NSs. 
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1. Introduction 

In recent years, great progress has been made in observa- 
tions of thermal radiation from middle-aged isolated NSs. 
Main observational results are summarized in Fig. |l| which 
shows error bars of the effective surface temperature, T s °° , 
as measured by a distant observer, versus stellar age t for 
eight sources. The youngest three are radio-quiet NSs in 
supernova remnants, the oldest is also a radio-quiet NS, 
while the others manifest themselves as radio pulsars. NS 
ages are mainly pulsar spindown ages or the estimated 
supernova ages. The age of the Vela pulsar is taken ac- 
cording to Lyne et al. (|1996 ), and the age of RX J 185635- 
3756 (lgi [yr] = 5.95) is taken according to Walter fl200l| ). 
The values of T s °° come fr om th e following sources: RX 
■10822- 43 — Zavlin et al. ( |l999[ ), IE 1207-52 — Zavlin 
et al. ([199SD , RX J0002+62 — Zav lin fc Pavlov ( |l999| ), 
PSR 0833-45 (Vela) — P avlov et al. ( |200lD , PSR 0656+14 

- Possenti et al fll996j ), PSR 0630+178 (Gemin ga) - 
Halpern & Wang ( |l997j ), PSR 10 55-52 — Ogelman ( |1995| ), 
and RX J1856-3754 — Walter (|200l|). For all the sources 



obtained from the observed X-ray spectra using hydro- 
gen atmosphere models. They are more consistent with 
other information on these sources (distances, hydrogen 
column densities, inferred neutron star radii, etc) than 
the blackbody model of NS radiation. On the contrary, 
the blackbody model is more consistent for older sources, 
and we present the blackbody values of T s °° for the older 
sources. We do not present observational data for some 
other sources. For instance, we do not cons ider a possible 
NS candidate in Cas A (e.g., Pavlov et al. 2000 ) because 
its nature is not clear. 

We will interpret the observational data with the mod- 
els of cooling NS. NSs of age t SS 10 6 yr are known to cool 
mainly via neutrino emission from their interiors, while 
the older NSs cool mainly via photon emission from the 
surface. We use our fully relativistic nonisothermal cooling 
code constructed by O. Gnedin and described in Gnedin 



but the Vela pulsar and RX J1856-375 4 thes e data are pre- 
sented in Table 3 in Yakovlev et al. ( |1999[ ). For the Vela 
pulsar T s °° = 0.68 ± 0.03 MK at lcr confidence level. For 
RX J1856-3754 we take T s °° = 0.57 ± 0.06 MK (the cen- 
tral value is mentioned by Walter; we supplement it with 
an estimated error for a joint fit of the optical plus X-ray 
data). The values of T s °° for the four youngest sources are 



et al. (2001). We employ the simplest NS models assuming 
that the NS cores consist of neutrons (n) , protons (p) and 
electrons. We adopt a moderately stiff equation of state 
(EOS) of this matter proposed by Prakash et al. (1988) 
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(their model I with the compression modulus of saturated 
nuclear matter K — 240 MeV). Other physics input is 
the same as described in Gnedin et al. (2001). In particu- 
lar, we place the core-crust interface at p = 1.5 x 10 14 
g cm -3 . The outer heat blanketing NS envelope is as- 
sumed to be made of iron (although the atmosphere may 
contain light elements) neglecting the effects of surface 
magnetic fields; the relation between T s and the internal 
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NS temperature is taken from Potekhin et al. (1997). The 
validity of such approach is discussed, e.g., by Yakovlev 
et al. (1999). The maximum mass of our NS models is 
M max = 1.977 M Q (with the radius of 10.754 km) and 
the maximum central density is p c ,max = 2.575 x 10 15 g 
cm -3 . For the given EOS, powerful direct Urea process of 



neutrino emission (Lattimer et al. 1991) is forbidden at 
M < M D = 1.358 M© (p c < p D = 7.851 x 10 14 g cm" 3 ). 
If so, a NS undergoes a slow cooling mainly via modified 
Urea process of neutrino emission. A NS with M > Md 
possesses a central kernel, where direct Urea process is 
open (p > pd); its cooling is fast. 

The values of T s °° for IE 1207-52, RX J0002+62, and 
PSR 0656+14 may be consistent with slow cooling of 
non-superfluid NSs. RX J0822-43, PSR 1055-52 and RX 
J1856-3754 are too hot for slow cooling of a non-superfluid 
NS with the heat-blanketing envelope made of iron (al- 
though RX J0822-43 is consistent with the bla nketin g en- 
velope made of light elements, Potekhin et al. 1997 ). The 
lower values of T s °°, especially for the Vela and Geminga 
pulsars, are intermediate between slow and fast cooling of 
non-superfluid NSs. In non-superfluid NSs, the transition 
from slow to fast cooling with increasing M occurs in a 
very narrow mass range 1 < M/ Mp ^ 1.003 (e.g., Page 
& Applegate 1992 , Yakovlev et al. 2001 ). In principle, one 
can explain the Geminga and Vela data by cooling of non- 
superfluid NSs but the probability that NS masses fall in 
this very narrow mass range is too low to accept such 
explanation. 

We will show that the observational data in Fig. [I] can 
be explained by the models of NSs with superfluid cores. 



2. Nucleon superfluidity 

We adopt traditional assumption that neutrons undergo 
singlet-state Cooper pairing in the NS crust and triplet- 
state pairing in the core, while protons undergo singlet- 
state pairing in the core. Microscopic calculations of the 
critical temperatures of the neutron and proton superflu- 
idities (SFs), T cn and T cp , depend sensitively on the model 
of nucleon-nucleon interaction and many-body theory em- 
ployed (see, e.g., Yakovlev et al. 1999 , for references). All 
models predict pronounced density dependence of T cn and 
T cp . In particular, T CD (p) for the singlet-state neutron SF 
has maximum at subnuclear densities in the crust and 
vanishes at p ~ 2 x 10 14 g cm -3 while T cn (p) for the 
triplet-state neutron SF grows up at subnuclear density, 
reaches maximum at p — (2 — 3) po (po — 2.8 x 10 14 g 



cm is the saturation density of nuclear matter) and de- 
creases with p, vanishing at p ~ (3 — 5) x 10 15 g cm' 3 . T cp 
also has maximum at several po and vanishes at higher 
p. At p ~ (2 — 3) po the proton SF is typically stronger 
than the neutron one. SF of n and p reduces the emis- 
sivity of neutrino reactions and affects heat capacity of n 
and p. Moreover, SF initiates a specific neutrino emission 
assoc iated with Cooper pairing of nucleons (Flowers et al. 
19761) . In this way nucleon SF becomes a strong regulator 



of NS cooling. The appropriate physics input is described 
in Yakovlev et al. fl200l| ). 

SF of neutrons in the crust affects cooling of young NSs 
(t 100 yr), in which the internal thermal relaxation is 
not over, and it does not affect cooling of middle-aged NSs 
we are interested in. To be specific, we neglect neutron SF 
in the crust. As for T cn (p) and T cp (p) in the NS core, we 
will not rely on any particular microscopic model but will 
try to find the models consistent with the observations. We 
parameterize the dependence of T c n on the nucleon (N=n 
or p) wavenumber k = /sen = (37t 2 ?in) 1 ^ 3 (measured in 
feu -1 ) as 
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for fc < k < fc 2 , and T c n = for k < fc and k > fc 2 , T , 
k ,. . .k 3 being free parameters. T regulates maximum of 
7cn(p), ko and k 2 determine low-density and high-density 
?cn(p) cutoffs, while k\ and k-$ govern the shape of T c n(p). 
Let us present also our simple fit of the density dependence 
of the nucleon number density tin (in fm~ 3 ) for the given 
EOS: 
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a Pi<i/( l + cpi4 + dpl i ), 
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where p u = p/(10 14 g cm' 3 ); a = 0.1675, b = 1.8185, 
c = 2.0288, d = 0.02444 for N=n, and a = 0.0006823, 
b = 2.6767, c = 0.1946, d = 0.01604 for N=p. 

At the first step we neglect neutron SF in the core and 
adopt T cp (p) displayed in Fig. ^. It is given by Eqs. ((TJ) 
and (|) with T = 2.029 x 10 10 K, k Q = 0, h = 1.117 
fin" 1 , fc 2 = 1.241 fin -1 , k 3 = 0.1473 fin" 1 . The maximum 
T cp = 6.78 x 10 9 K takes place at p = 5.74 x 10 14 g cm" 3 
and T cp vanishes at p = 9.55 x 10 14 g cm' 3 . This proton 
SF is constructed artificially but is typical for SFs pro- 
vided by microscopic calculations. In a not too massive 
NS (M < 1.45 M , T cp (p c ) ^ 10 9 K) it will appear at 
the initial stage of NS cooling (t ^ 100 yr). Accordingly, 
the peak of neutrino emission due to Cooper pairing of p 
will take place at this early stage making almost no affect 
on the cooling curves (solid lines). Moreover, the internal 
temperature of the middle-aged NS will be much lower 
than T cp . The proton SF will switch off the proton heat 
capacity and it will strongly suppress neutrino reactions 
involving protons, particularly, modified and direct Urea 
processes. However, the total heat capacity is mainly pro- 
vided by neutrons being almost unaffected by proton SF. 

The upper cooling curve is for the 1.3 M Q NS (with 
the radius of 13.04 km) but it would be exactly the same 
for all M from Md down to ~ 1.1 M Q (similar insensitiv- 
ity to M is discussed, e.g., in Yakovlev et al. |1999| ). For 
t ~ (10 2 — 10 5 ) yr this cooling curve goes noticeably higher 
than in a non-superfluid NS since SF suppresses neutrino 
emission and slows cooling of such a NS. This allows us 
to explain high surface temperature of RX J0822-43 by 
cooling of a superfluid NS with the heat blanketing enve- 
lope made of iron. Similar to the nonsuperfluid NS mod- 
els, we have faster cooling at M > Md- However, the SF 
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Fig. 1. Observational data (error bars) on surface temper- 
atures of eight NSs (see the text) as compared with theo- 
retical cooling curves obtained for proton and neutron SFs 
from Fig. g. All cooling curves (except dot-and-dashed 
ones) are calculated for the same proton SF. Solid lines 
- no neutron SF for NS models with masses (from top to 
bottom) 1.3, 1.36, 1.4, 1.46, 1.47, and 1.5 M . Dashed lines 
and dotted lines correspond to neutron SFs for M = 1.3 
and 1.4 M . Dashed-and-dot lines are for non-superfluid 
1.3 M Q and 1.5 M Q NSs. 



suppression of neutrino emission strongly smoothes tran- 
sition from slow to fast cooling and allows us to explain 
the low temperatures T s °° of the Vela and Geminga pul- 
sars by cooling of NSs with M w 1.47 M . With further 
increase of M to M max we obtain almost the same fast- 
cooling curves as in the absence of SF. These results are 
not sensitive to the specific choice of T cp (p) near the max- 
imum (at p ^ pu) as long as T cp (p) is high (say, ^3x 10 9 
K) for p <, pY). However, the transition from slow to fast 
cooling with increasing M is very sensitive to the decreas- 
ing slope of T cn (p) at higher p. By shifting this slope to 
higher p we would obtain higher masses of the Vela and 
Geminga pulsars. 

Our next step is to discuss the effects of neutron SF. 
Following the results of microscopic calculations we as- 
sume that this SF is generally weaker than the proton 
one in the NS core. If T cn <, 3.5 x 10 8 K, this SF ap- 
pears in the NS core at t > 100 yr (for M <; 1.45 M Q ). 
The NS core is then isothermal, and the temperature de- 
creases uniformly over the core in the course of cooling. 
Thus the SF appears first in a spherical layer at those val- 
ues of p = pep which correspond to maximum of T cn (p), 
Tcn(pcp) = T™ aK . Before that time the cooling curves 
are exactly the same (solid lines in Fig. |l|) as in the ab- 
sence of neutron SF. Once the neutron SF is switched on, 
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Fig. 2. Density dependence for one model of the proton 
critical temperature (solid line) and two models of neutron 
critical temperatures (dots and dashes) used in cooling 
simulations (Fig. Q). Vertical dot-and-dashed lines show 
the central densities of 1.1, 1.3, 1.4, 1.5 and 1.6 Mq NSs. 



powerful neutrino emission due to Cooper pairing of n 
starts to operate. If pep is noticeably lower than the NS 
central density one can calculate the neutrino luminosity 
Lqp due to this Cooper pairing semianalytically. Its tem- 
perature dependence can be written as Lqp{T) = Lq 1{t), 
where Lq is a normalization constant (which depends on 
T C n(p) and NS model) while 1(t) is a universal function of 
r = T/T™ ax (with I oc (1 - r) 3 / 2 at (1 - t) < 1; I oc t 8 
at r < 1, and with maximum at t w 0.795). This neu- 
trino emission cools the star very quickly producing sharp 
breaks of the cooling curves. The breaks are especially 
pronounced if the neutron SF appears at t ~ 10 4 - 10 6 
yr. Soon after the break onsets the cooling curves become 
even steeper due to reduction of the neutron heat capac- 
ity at T <C Ten and due to transition to photon cooling 
stage. The effect is illustrated in Fig. [I] which shows cool- 
ing curves (dots and dashes) for two models of neutron SF 
(dots and dashes) displayed in Fig. ^ and parameterized 
byEqs. (Q) and (§). In both cases ko = 1 fm 1 . The dotted 
line in Fig. |^ corresponds to the following set of parame- 
ters: T = 2.915 x 10 8 K, ki = 0.5664 fm" 1 , k 2 = 2.767 
fm -1 , &3 = 0.2965 fm -1 . Accordingly, T cn has maximum 
1.93 x 10 8 K at p = 5.62 x 10 14 g cm' 3 and vanishes at 
p = 1.80 x 10 15 g cm~ 3 . The dashed line corresponds to 
To = 6.461 x 10 9 K, ki = 1.961 fin -1 , fc 2 = 2.755 far 1 , 
fa = 1.30 fm" 1 . Thus, T cn has maximum 3.41 x 10 8 K at 
4.52 x 10 14 g cm" 3 and vanishes at p = 1.78 x 10 15 g 

3 



S3 
P = 

cm 



The dashed-curve SF has somewhat larger maximum 
T cn and appears at earlier cooling stage. As seen from Fig. 
0, any mildly strong neutron SF destroys our theoretical 
explanation of observational data producing pronounced 
breaks of the cooling curves. It is evident that the weaker 
SF of neutrons violates theoretical interpretation of older 



4 



A.D. Kaminker et al.: Nucleon superfluidity versus masses in neutron stars 



NSs. It would be impossible to explain rather high ob- 
served surface temperature of old PSR 1055-52 and RX 
J 1856-3754 by our NS models with sufficiently strong SF 
of neutrons (regardless of SF of protons). The NS heat ca- 
pacity would be strongly reduced, neutrino emission due 
to Cooper pairing of neutrons would be important, and 
the NS would be much colder. Our interpretation of all 
sources in Fig. ^ would not be violated only by a weak 
neutron SF with maximum T cn < 10 8 K. This result is in 
line with low values of T cn obtained earlier (Yakovlev et 
al. 1999) using simplified models of T cn consta nt thr ough- 
out the NS cores. Note that Yakovlev et al. (11999D were 
unable to explain rather high values of T s °° for PSR 1055- 
52. We can provide such an explanation using NS models 
with somewhat different EOS and more advanced physics 
input in the NS crust. 



3. Summary 

We have proposed an interpretation of the observations us- 
ing a simple model of cooling NSs. The model is consistent 
with observational data if the proton SF in dense matter 
is rather strong, with the maximum of T cp <; 5 x 10 9 K 
at p ~ (2 — 3) po but decreases sharply at p <^ po- On the 
other hand, the neutron SF must be weak, with the max- 
imum T cn < 10 8 K. Fixing the density profiles of T cp (p) 
and T cn (p) we can explain the observational data by cool- 
ing of NSs with the same EOS and SF properties but with 
different masses. Our interpretation does not require any 
reheating mechanism in cooling isolated NSs (see, e.g., 
Yakovlev et al. 1999, for references). 

A strong SF of protons is required to smooth transi- 
tion from slow to fast cooling with increasing NS mass. 
Otherwise it would be difficult to explain rather low sur- 
face temperatures of the Vela and Geminga pulsars. A 
weakness of neutron SF is necessary to avoid sharp falls of 
the NS surface temperatures after onset of neutrino emis- 
sion due to Cooper pairing of neutrons. Weak neutron SF 
is in favor of a not too soft EOS in the NS core (softness 
would mean strong attractive nn interaction and therefore 
strong neutron pairing). Strong proton SF is in favor of 
not too large symmetry energy at high densities (too large 
symmetry energy would mean large proton number den- 
sity which would suppress proton pairing). On the other 
hand, the symmetry energy must be not too small to open 
direct Urea process at p > po • 

The adopted density dependence T cp (p) is not unique. 
We could choose different T cp (p) profiles satisfying the 
above criteria and interpret observations in the same man- 
ner by cooling of NSs with somewhat different masses. 
Although we have used one particular EOS in the NS core, 
about the same results would evidently be obtained for a 
variety of EOSs. 

In our model, the colder NSs, Vela and Geminga, 
are more massive (M ~ 1.45 M©), while the hottest RX 
J0822-43 and warm and old PSR 1055-52 and RX J1856- 
3754 are the least massive (M ~ (1.1 -1.3) M Q ) of all NSs 
included in the analysis. The mass range obtained is in 



good agreement with the well-known range of masses of 
radio pulsars in binary systems (Thorsett & Chakrabarty 
1999). Note that in our model neutrons in the NS cores 
remain nonsuperfluid for t 10 6 yr. Furthermore, mas- 
sive NSs, M J; 1.5 M Q , must undergo fast cooling and be 
cold, T s °° <, 3 x 10 5 K (the lowest solid curve in Fig. Q). 
Thus, we could explain the existence of cold middle-aged 
NSs, if such objects were discovered, as massive NSs. On 
the other hand, it would be easy to modify the model to 
avoid prediction of these cold NSs. This can be done by 
choosing the EOS in the NS core in which p® is close to 
Pc,max- Another possibility is to adopt the T cp (p) profile 
with T cp (p Cjmax ) £ (2 - 3) x 10 9 K. 

Our model may seem oversimplified because it ne- 
glects possible presence of other particles (muons, hyper- 
ons, quarks) in the NS cores. However, despite its simplic- 
ity, it can explain observations of isolated NSs with only 
one adjustable parameter - NS mass. We are planning to 
analyze more complicated models in future publications. 
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